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■ ABSTRACT. Given a self-adjoint involution 7 on a Hilbert space f), we consider ay-self-adjoint 

— ■ operator L = A + V on $j where A is a possibly unbounded self-adjoint operator commuting 

with J and V a bounded ./-self-adjoint operator anti-commuting with J. We establish optimal 
ryy . estimates on the position of the spectrum of L with respect to the spectrum of A and we obtain 

norm bounds on the operator angles between maximal uniformly definite reducing subspaces of 
the unperturbed operator A and the perturbed operator L. All the bounds are given in terms of the 
norm of V and the distances between pairs of disjoint spectral sets associated with the operator L 
0^ ' and/or the operator A. As an example, the quantum harmonic oscillator under a ^-symmetric 

ry) , perturbation is discussed. The sharp norm bounds obtained for the operator angles generalize 

• ■ the celebrated Davis-Kahan trigonometric theorems to the case of /-self-adjoint perturbations. 

-I— > 

S: 

1. INTRODUCTION 

Let f) be a Hilbert space and / a self-adjoint involution on 9), that is, J* =J and J 2 = I, where 
J ^ I and / denotes the identity operator. A linear operator L on S) is called /-self-adjoint if the 
product JL is a self-adjoint operator on $), that is, (JL)* = JL. 

In this paper we consider a /-self-adjoint operator L of the form L = A + V where A is a 
(possibly unbounded) self-adjoint operator on Sj commuting with / and V a bounded /-self- 
adjoint operator anti-commuting with J. Since the involution J is both unitary and self-adjoint, 
qs, ' its spectrum consists of the two points + 1 and — 1 and hence 



o 



/=E 7 ({+1})-E 7 ({-1}), (1.1) 



where E/({±1}) denote the corresponding spectral projections of /. Thus, the involution J 
^ | induces a natural decomposition of the Hilbert space Sj into the sum 

£ = -fto©-fti (1.2) 

of the complementary orthogonal subspaces 

^ = RanE / ({+l}), tfi =RanE 7 ({-l}). (1.3) 

Our assumptions on the operators A and V imply that they are nothing but the diagonal and 
off-diagonal parts of L with respect to the decomposition (11.21 ): 

^° ^ Y Dom(A)=Dom(A )©Dom(Ai)c£o©£i, (1.4) 
B Q V Dom(V)=j3o©£i; (1-5) 
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here the entries Aq =A\^ and A\ = A| . are self-adjoint operators on S)q and respectively, 
and B = V |^ G &($)i,Sjo) is bounded. Thus the operator L may be viewed as an off-diagonal 
bounded /-self-adjoint perturbation of the block diagonal self-adjoint operator matrix A. 

A powerful tool to study operators L admitting a block operator matrix representation with re- 
spect to a self-adjoint involution J is furnished by indefinite inner product spaces. This idea was 
first used in ||26l to prove a general theorem on block-diagonalizability for a ./-accretive operator 
A and a self-adjoint perturbation V, with application to Dirac operators. The main ingredient 
of this approach is to show that the perturbed reducing subspaces are maximal uniformly posi- 
tive and negative, respectively, with respect to the indefinite inner product. As a consequence, 
they admit graph representations by angular operators which measure the deviation between the 
unperturbed and the perturbed reducing subspaces. 

In the situation considered in the present paper, the self-adjoint involution / induces an in- 
definite inner product by means of the formula 

\x,y] = (Jx,y), x,ye&. (1.6) 

The Hilbert space Sj equipped with the indefinite inner product (11.61 ) is a Krein space which we 
denote by 8, assuming that R stands for the pair {$),J}. Note that every /-self-adjoint operator 
on fj is a self-adjoint operator on the Krein space in particular, the operators A, V, and 
L = A + V are self-adjoint operators on In the Krein space a (closed) subspace £ C 8. is 
said to be uniformly positive if there exists a y > such that 

[x,x] > 7 ||^|| 2 for every x G J?, x / 0, (1.7) 

where || • || denotes the norm on $). The subspace £ is called maximal uniformly positive if it is 
not a proper subset of another uniformly positive subspace of 8. Uniformly negative and max- 
imal uniformly negative subspaces of & are defined in a similar way, replacing the inequality in 
(11.71) by [x,x] < — 7 ||x|| 2 . Direct sums of subspaces of & (or ft) that are /-orthogonal (i.e. ortho- 
gonal with respect to the inner product [•,•]) are denoted with "[+]". Further definitions related 
to Krein spaces and linear operators therein may be found, e.g., in ll23ll . iPTTl . lPT6ll . or Q • 

The subspaces f)o and Sji, which simultaneously reduce A and J, are maximal uniformly 
positive and maximal uniformly negative, respectively, with respect to the inner product (11.61 ) 
induced by J. Throughout this paper, we assume that also the perturbed operator L = A + V 
possesses a maximal uniformly positive invariant subspace S)' . Then the complementary J- 
orthogonal subspace f)\ = is invariant for L as well; hence both f>Q and Sj[ are automat- 
ically reducing subspaces for L and the spectrum of L is purely real (see, e.g., Corollary 12. 121 
below). 

The main goal of this paper is to establish bounds on the position of the reducing subspaces 
Sj'q or $j\ of the perturbed operator L = A + V relative to the subspaces Sjqov Sji. The bounds are 
given in terms of the norm of the perturbation V and of the distances between the unperturbed 
spectra 

a = spec(A| flo ) and o { = spec(A|^) (1.8) 

of A and/or the perturbed spectra 

do = spec (L I g/ ) and a[ = spec (L| fl , ) (1.9) 

of L in their respective maximal uniformly definite reducing subspaces. 

We describe the mutual geometry of the maximal uniformly definite reducing subspaces of 
the unperturbed and perturbed operators A and L = A + V by using the concept of operator angles 
between two subspaces of a Hilbert space (for a discussion of this concept and references see, 
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e.g., HH). Recall that the operator angle 0(f),-,f)-) between Sji and measured relative to fit 
is given by (see, e.g., Ifl9l0 



0, = ©(«,, j5{) = arcsm^/ A -P j 5 t P j5 /| j5( , f = 0,1, (1.10) 

where 1%, denotes the identity operator on fj,, and and stand for the orthogonal pro- 
jections in fj onto fj,- and f)-, respectively. By definition, the operator angle &($)(,$)';) is a 
non-negative operator on fj; and 

||0(^,^)|| =maxspec(0(^,^)) < n/2. 

The main tool we use to estimate the operator angles 0,- = 0(f},-,f}-) is their relation to 
solutions of the operator Riccati equation 

KAq —AyK + KBK = —B* , (1.11) 



where the coefficients Aq, A\, and B are the entries of the block matrix representations (11.41) and 
(11.51 ) of the operators A and V. In fact, given a maximal uniformly positive reducing subspace 
#q of L = A + V, there exists a unique uniformly contractive solution K (\\K\\ < 1) to the 
Riccati equation (11.111 ) such that fj is the graph of K; the maximal uniformly negative reducing 
subspace S)\ of L, which is /-orthogonal to $j' , is the graph of the adjoint of K. Since \\K\\ < 1 
and \K\ = tan0(fj,,.f)-) (see Remark 12761 and Lemma 12781 below), the operator angle always 
satisfies the two-sided inequality 

O<0,<7T/4, i = 0,1. (1.12) 

By establishing tighter norm bounds on the uniformly contractive solution K of (11.111 ). we thus 
obtain tighter norm bounds for the operator angles (1 1 - 10b . 

Sufficient conditions guaranteeing the existence of maximal uniformly definite reducing sub- 
spaces for the operator L = A + V may be found, e.g., in ||6l and |[34ll35l . In particular, one of 
the main results of O is as follows. Here and in the sequel, by conv(a) we denote the convex 
hull of a Borel set a C R. 

Theorem 1.1 (|61, Theorem 5.8 (ii)). Assume that the spectral sets Co and 0\ are disjoint, i.e. 

d :=dist(ob,oi) > 0, (1.13) 
that one of these sets lies in a finite or infinite gap of the other one, i.e. 

conv(a,-) n G\-i = or <7 ! -nconv(cTi_ 1 -) =0 for some i = 0,1, 
and that \\V\\ <d/2. Then 

spec(L) = OqG<j[, 

where the (disjoint) sets Cq C R and ffjcK lie in the closed \\V\\/2-neighbourhoods of the sets 
Ob and G\, respectively. The complementary spectral subspaces $)'q and f)\ of L associated with 
the spectral sets Cq and o[ are maximal uniformly positive and maximal uniformly negative, 
respectively, and satisfy the sharp norm bound 

(\ 2||V||\ 
tan©,- < tanh - arctanh , i = 0, 1, (1.14) 



or, equivalently, 

|| sin 20; || < ^P-, i = 0,l. (1.15) 
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The bound (11.14b relies on the disjointness of the spectral sets Ob and Oy of the unperturbed 
operator A and involves the distance between Oq and Oi . Therefore, this bound (as well as the 
other bounds from [6, Theorem 5.8]) is an a priori estimate. In the present paper, we establish 
bounds on the operator angles 0, that involve at least one of the perturbed spectral sets <7q and 
a[. In general, for these new bounds to hold, the disjointness ( 11.13b of the sets Go and G\ is not 
required at all. 

Our first main result is a semi-a posteriori bound on the operator angles 0, involving the dis- 
tances dist((7o, o[) and/or dist(<7i,<7o) between one unperturbed and one perturbed spectral set. 

Theorem 1.2. Suppose that L has a maximal uniformly positive reducing subspace fy' Q in the 
Krein space .ft = {?),J}. 

(i) If for some i = 0, 1 the spectral sets a, and a[_ ; - are disjoint, i.e. 

5 i :=dist(a i -,a 1 / _ ; ) > 0, (1.16) 
then the operator angles &j satisfy the bound 

tan©; < forboth j = 0,1. (1.17) 

2 Oi 

(ii) If in addition, one of the sets (7, and o[_ i satisfying (11.16b lies in a finite or infinite gap 
of the other one, i.e. 

conv(Gi) n o[_j = or a, nconv(a 1 '_ ! ) = 0, (1-18) 

then we have the stronger estimate 

||tan0y|| < forboth j = 0,1. (1.19) 

Our second main result is a completely a posteriori estimate since it only involves the distance 
between the spectral sets Oq and a[ associated with the perturbed operator L. 

Theorem 1.3. Suppose that L has a maximal uniformly positive reducing subspace f)' Q in the 
Krein space ^ = {fj,/}. 

(i) If the spectral sets o' Q and o[ are disjoint, i.e. 

8 :=dist(a o ,a0 > 0, (1.20) 
then the operator angles &j satisfy the estimate 

I|tan0/|| < -S for both j = 0,1. (1.21) 
2 § 

(ii) If, in addition, for some i = 0, 1 the set o[ is bounded and lies in a finite or infinite gap 
ofa[_ v i.e. 

conv(a/)na 1 '_,. = 0, (1.22) 
then we have the stronger estimate 

WW 

tan 0; 1| < ; " " forboth j = 0,1. (1.23) 

<5 2 + ||V|| 2 
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(iii) Furthermore, if both spectral sets o' {) and o[ are bounded and subordinated, i.e. 



conv(ao) n conv(oi) = 0, 



(1.24) 



then we have the even stronger estimate 



tan 2®; 1 1 < 



2\\V 



for both 7 = 0,1. 



(1.25) 
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The bounds d 1 - 19b and (11.251 ) as well as the bound (11.231 ) in the case of a finite gap are 
optimal (see Remarks I4.lj - I4.3l ). Moreover, the sharp a priori bound (11.141) turns out to be a 
corollary either to Theorem ll.2l (ii) or to Theorem ll.3l (iii) (see Theorem 16.41 and Remark 16.81 
respectively). 

The semi-a posteriori bounds of Theorem 11.21 and the completely a posteriori ones of The- 
orem 11.31 complement the a priori norm bounds on the variation of spectral subspaces for J- 
self-adjoint operators proved in |6, Theorem 5.8]. The sharp norm bounds of these theorems 
represent analogues of the celebrated trigonometric estimates for self-adjoint operators known 
as Davis-Kahan sin®, sin2®, tan®, and tan2® theorems (see [13] and the subsequent papers 
0[]j3|2Ol|2l][3Tl): ^ e bound (11.151 ) may be called the a priori sin2® theorem for /-self-adjoint 
operators; the bounds d 1 - 19b and (11.231 ) may be called the semi-a posteriori and completely a 
posteriori tan® theorems, respectively; the bound (11.251 ) may be called the a posteriori tan 2® 
theorem. 

The plan of the paper is as follows. In Section |2] we give necessary definitions and recall 
some basic results on the block diagonalization of /-self-adjoint 2x2 block operator matrices. 
In Section [3] we establish several semi-a posteriori and completely a posteriori norm bounds on 
uniformly contractive solutions to operator Riccati equations of the form (11-1 lb - Using these 
results, we prove both Theorems 11.21 and 11.31 in Section 01 Assuming that the spectral sets 
(11.81) do not intersect and ||V|| < ^dist(ao,ai), in Section [5] we obtain sharp estimates on the 
position of the isolated components of the spectrum of L = A + V confined in the closed ||V||- 
neighbourhoods of the sets Oq and Oi . In this section, we also establish bounds on the spectrum 
for more general 2x2 block operator matrices that need not be /-self-adjoint. In Section [6j 
we combine Theorems 11.21 and 11.31 with the spectral estimates of Section [5] and discuss the 
emerging a priori norm bounds on variation of the spectral subspaces of a self-adjoint operator 
on a Hilbert space under /-self-adjoint perturbations. Finally, in Section[7]we apply some of the 
bounds obtained to the Schrodinger operator describing an A^-dimensional isotropic harmonic 
oscillator under a ^^-symmetric perturbation (see e.g. (H); here the parity operator & plays 
the role of the self-adjoint involution J (see |[2ll6l l27l . 

The following notations are used throughout the paper. By a subspace of a Hilbert space we 
always mean a closed linear subset. The identity operator on a subspace (or on the whole Hilbert 
space) 9JT is denoted by Ajjt; if no confusion arises, the index 9JT is often omitted. The Banach 
space of bounded linear operators from a Hilbert space fj to a Hilbert space $)' is denoted by 
and by SS{¥$) if S) = Sj'. The symbol U is used for the union of two disjoint sets. By 
O r (L), r > 0, we denote the closed r-neighbourhood of a Borel set £ in the complex plane C, 
i.e. O r (L) = {z G C dist(z,E) < r}. By a finite gap of a closed Borel set C K, a / 0, we 
understand an open interval (a,b), — °° < a < b < oo, such that on(a,b) =0 and a,b £ a; an 
infinite gap of a is a semi-infinite interval (a, b) such that a n(a,b) = and either a = — oo, 
\b\ < oo, and b G a or \a\ < °°, a £ a, and b = °°. 
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2. Preliminaries 

In this section we recall some results on the block diagonalization of /-self-adjoint operator 
matrices in terms of solutions to the related operator Riccati equations and on norm bounds for 
solutions to operator Sylvester equations. We also recall a couple of statements on maximal 
uniformly definite subspaces of a Krein space. For notational setup we adopt the following 

Assumption 2.1. Let J be a self-adjoint involution on a Hilbert space ft, J ^ I, and let fto 
and ft\ be the spectral subspaces (11.31 ) of J. Also assume that A is a (possibly unbounded) 
self-adjoint operator on ft diagonal with respect the decomposition (11.21 ), which means that fto 
and ft\ are the reducing subspaces of A and the representation (11.41 ) holds withAo and Ai the 
self-adjoint operators on fto and ft h respectively. Let V be a bounded operator on ft admitting, 
relative to (11.21 ). the representation (11.51 ) where B G 3§(fti,fto). Finally, let 

L = A + V=(_2* *Y Dom(L) = Dom(A). (2.1) 

With a block operator matrix L of the form (12.11 ) we associate the operator Riccati equation 
(11.111) where K is a linear operator from fto to ft \ . There are different concepts of solutions to 
such an equation; here we recall the notion of weak and strong solutions (see (3l 0]). 

Definition 2.2. Assume that Assumption 12. II is satisfied. A bounded operator K G £$(fto,fti) 
is said to be a weak solution to the Riccati equation (11.111) if 

(KA Q x,y) - (Kx,A\y) + (KBKx,y) = -(B*x,y) for all x G Dom(A ), y G Dom(A* 1 ). 

A bounded operator K G &(fto,ft\) is called a strong solution to the Riccati equation (11.111 ) if 

Ran(£| Dom(Ao) ) cDom(Ai) (2.2) 

and 

KA x -AiKx + KBKx = -B*x for all x G Dom(A ) . (2.3) 

Remark 2.3. Obviously, every strong solution K G 38 {fto, ft \) to the Riccati equation (11.111 ) is 
also a weak solution. In fact, the two notions are equivalent by Lemma 5.2]: every weak 
solution K G 38 {fto, ft \) of the operator Riccati equation (11.111) is also a strong solution. 

Remark 2.4. With the block operator matrix (12.11) . one can also associate the operator Riccati 
equation 

K'Ai -AqK' - K'B*K' = B (2.4) 
where K' is a linear operator from f)i to fjo- From Definition 12.21 it immediately follows that 
K G &($)o,Sj\) is a weak (and hence strong) solution to (11.111) if and only if K' = K* is a weak 
(and hence strong) solution to (12.4I ). 

Definition 2.5. Let 9Jt be a subspace of the Hilbert space ij, SOT- 1 = S) Q 9JT its orthogonal 
complement, and K a bounded linear operator from 9Jt to 9JI . Denote by P^t and P w ±_ the 
orthogonal projections in S) onto the subspaces 971 and 3Jt , respectively. The set 

<S{K) = {x G ft | P m LX = KP m x} 

is called the graph subspace associated with the operator K. 

Remark 2.6. If a subspace © C ft is a graph (5 = &(K) of a bounded linear operator K G 
=5^(271,971 ), then K is called the angular operator for the (ordered) pair of subspaces 97T and 
©; the usage of this term is explained by the equality (see lfT8l : cf. lTT3l and lTT4l ) 

\K\ =tan0(97t,0), (2.5) 
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where \K\ is the modulus of K, \K\ = \/ r K*K, and &(VJi, (3) is the operator angle between the 
subspaces 9JT and <3 measured relative to the subspace 9Jt (see definition (11.101) ). 

It is well known that strong solutions to the Riccati equations (11.111 ) and (12.41 ) determine 
invariant subspaces for the operator matrix L by means of their graph subspaces (see, e.g., [3] 
and 11241 ) . More precisely, the following correspondences hold (see, e.g., (6l Lemma 2.4]). 

Lemma 2.7. Assume that Assumption I2.il holds. Then the graph @(K) of an operator K G 
3£l($jo,$)l) satisfying (12.21 ) is an invariant subspace for the operator matrix L if and only if K 
is a strong solution to the operator Riccati equation (ll.lll ). Similarly, the graph < £(K') of an 
operator K' G ,i^o) is an invariant subspace for L if and only ifK' is a strong solution to 

the Riccati equation (12.4I ). 



The next two statements are well-known facts in the theory of spaces with indefinite metric 
(see, e.g., [7J Section 1.8, in particular, Corollaries 1.8.13 and 1.8.14]). 

Lemma 2.8. A subspace £ is a maximal uniformly positive subspace of the Krein space R if 
and only if there is a uniform contraction K G &($)o,f)i) (i.e. \\K\\ < I) such that £ is the 
graph @(K) of the contraction K. Similarly, a subspace £' is a maximal uniformly negative 
subspace of the Krein space R if and only if £' is the graph ^(K 1 ) of a uniform contraction 
K' G &(f)iM. 

Lemma 2.9. Let £ be a maximal uniformly positive subspace of the Krein space Then the 
orthogonal complement £W f £ i n & i s a maximal uniformly negative subspace. If £ is a 
graph subspace, £ = &(K) with K G &($)o,$)i)> then £W is the graph of the adjoint of K, i.e. 
£W = &(K*), and £[+]£W = SL 



Many more details on Krein spaces and linear operators therein may be found in [22], [23], 
(LU, Hil or 0. 

The following sufficient condition for a /-self-adjoint block operator matrix of the form (12.11 ) 
to be similar to a self-adjoint operator on S) was proved in |6j; for the particular case where 
the spectra of the entries Aq and A\ are subordinated, say max spec (Aq) < minspec(Ai), closely 
related results may be found in d Theorem 4.1] and OOl Theorem 3.2]. 



Theorem 2.10 ([6], Theorem 5.2). Assume that L = A + V satisfies Assumption I2.il Suppose 
that the Riccati equation (11.1 lb has a weak (and hence strong) solution K G !r){) such 

that \\K\\ < 1. Then: 

(i) The operator matrix L has purely real spectrum and it is similar to a self-adjoint ope- 
rator on Sj. In particular, the equality 

L = TAT 1 (2.6) 

holds, where T is a bounded and boundedly invertible operator on f) given by 

I K* \ ( I K*K \ — ^/^ 

K I )\ I-KK* ) (2-7) 
and A is a block diagonal self-adjoint operator on f), 

A = diag(A ,Ai), Dom(A) = Dom(A ) ©Dom(Ai), (2.8) 
whose entries 

A = (1-K*K) { I 2 (A Q + BK)(1-K*K)-'I\ 

Dom(A )=Ran(i-^)V2| Dom(Ao) , ^ 
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and 

Ai = (I-KK*y/ z (Ai -B*K*)(l-KK*)- l l 2 , 

Dom(A 1 )=Ran(/-^)V2| Dom(Ai) , ( 2 -10) 

are self-adjoint operators on the corresponding Hilbert space components S)o and Sj\, 
respectively. 

(ii) The graph subspaces S)' = ^{K) and S)[ = @(K*) are invariant under L, mutually 
orthogonal with respect to the indefinite inner product (11.61 ). and 

The subspace S)' is maximal uniformly positive, while Sj\ is maximal uniformly nega- 
tive. The restrictions of L onto f)' and Sj\ are R-unitary equivalent to the self-adjoint 
operators Aq and A\, respectively. 

Remark 2.11. The requirement \\K\\ < 1 is sharp in the sense that if there is no uniformly 
contractive solution to the Riccati equation (ll.lll ). then the operator matrix L need not be similar 
to a self-adjoint operator at all; this can be seen, e.g., from (6] Example 5.5]. 

An elementary consequence of Theorem l2. lOl is the following property of maximal uniformly 
definite subspaces of /-self-adjoint operators L = A + V with self-adjoint A and bounded V. 

Corollary 2.12. Assume that L = A + V satisfies Assumption \2.1\ Suppose that L has a maximal 
uniformly positive (resp. negative) invariant subspace &o of R = {S),J}. Then &\ = 8^ is 
also an invariant subspace of L, which is maximal uniformly negative {resp. positive); the 
restrictions of L to &q and &\ are ^.-unitary equivalent to self-adjoint operators on the Hilbert 
spaces S)q and S)\, respectively. 

Proof. We give the proof for the case where is a maximal uniformly positive subspace; the 
proof for maximal uniformly negative $o is analogous. 

By Lemma 12781 £0 is the graph of a uniform contraction K : S)o — > £)i- Since ^0 is invariant 
under L, Lemma 12.71 shows that K is a uniformly contractive strong solution to the Riccati 
equation (11.1 lb . Now all claims follow immediately from Theorem l2.10l □ 

Riccati equations are closely related to operator Sylvester equations (also called Kato-Rosen- 
blum equations). In this paper we use the following well known result on sharp norm bounds 
for strong solutions to operator Sylvester equations (cf. 10 Theorem 4.9]). 

Theorem 2.13. Let Aq and Ai be (possibly unbounded) self-adjoint operators on the Hilbert 
spaces Sjq and respectively, and Y G 3$($Jo,$)i). If the spectra spec(Ao) and spec(Ai) are 
disjoint, i.e. 

d := dist(spec(Ao),spec(Ai)) > 0, 
then the operator Sylvester equation 

XA -A Y X = Y 

has a unique strong solution X S &(Sjo,Sji); the solution X satisfies the norm bound 

1 1 „ n\\Y\\ 

11*1 < 5 Ll; a") 

if in addition, one of the sets spec(Ao) and spec (A 1) lies in a finite or infinite gap of the other 
one, then X satisfies the stronger norm bound 

||X|| < 1H. (2.12) 
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Remark 2.14. The fact that the constant n/2 in the estimate ( 12.111 ) for the generic disposition of 
the sets spec(Ao) and spec (A i) is best possible is due to R. McEachin |[29l . The existence of the 
bound (12.121 ) for the particular case where one of the sets spec(Ao) and spec (A i) lies in a finite 
or infinite gap of the other one may be traced back to E. Heinz lPT5ll (also see O Theorem 3.2] 
and (6l Theorem 3.4]). For more details and references we refer the reader to (3l Remark 2.8] 
and Remark 4.10]. 

In the proofs of several statements below we will use the following elementary result, the 
proof of which is left to the reader. 

Lemma 2.15. Let <p be a scalar analytic function of a complex variable z whose Taylor series 

/ x ^ k 1 d k q>(0) 

<p(z) = £a^, % = ^-fF> * = 1,2,..., (2.13) 

k=o z 

is absolutely convergent on the open disc {z £ C : \z\ < r} for some r > 0. Let M G &($)i,fio) 
and N G £%($jo,$)i) be bounded operators with \\MN\\ < r and \\NM\\ < r. Then 

M(p(NM) = (p(MN)M, (2.14) 

where for a bounded linear operator T on a Hilbert space T with \\T\\ < r the value of (p(T) is 
defined by the series 

DO 

(p(T) = J £a k T k . 

k=0 

We also need the following auxiliary statement. 

Lemma 2.16. Assume that Assumption \2J] holds and suppose that the Riccati equation (11.111 ) 
has a weak (and hence strong) solution K G £$(f}i,$}o) such that \\K\\ < 1. Then 

Ran (*l D o m (A )) cDom ( A i) ( 2 -15) 

and 

KA y-A l Ky=-(I-KK*) l / 2 B*(I-K*K) l / 2 y for all y G Dom(A ), (2.16) 

where Aq and A{ are the self-adjoint operators given by ( 12.91 ) and (12.101 ). respectively. 

Proof. It is straightforward to verify that if K is a strong solution to the Riccati equation (11.111) . 
then 

KZ Q x - Z x Kx = -B*(I - K*K)x for all x G Dom(A ). (2.17) 
Since K is assumed to be a uniform contraction, Theorem l2.10l (i) applies and yields 

K(I - K*K)-^ 2 A {I - K*K) l l 2 x -(I- KK*)- l l 2 A x (I - KK*Y' 2 Kx 

= -B*(l-K*K)x for all x G Dom(A ). (2.18) 

Applying (/ - KK*) l l 2 from the left to both sides of (12481) and choosing x = (I- K*K)- l / 2 y 
with y G Dom(Ao), we arrive at the Sylvester equation 

XA a y - AiXy = Yy for all y G Dom(A ), (2.19) 

where 

X= (I-KK*) l/2 K(I-K*K)- l/2 , (2.20) 

Y = - (I-KK*) 1/2 B*(I-K*K) l/2 . (2.21) 
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Note that we have Ran(7 -^*^) _1/2 | Dom(Ao) = Dom(A ) by (1X91 ). Ran(£| Dom{Ao) ) C Dom(Ai) 
by d2!2t . and thus, by (12T01 ), 

Ran (^! D o m (Ao)) CDOm ( A l)- ( 2 - 22 ) 

Hence X is a strong solution to the Sylvester equation ( I2.19I ). 

Furthermore, the Taylor series (12.131 ) of the function <p(z) = (1 — z) 1 / 2 is absolutely con- 
vergent on the disc {z G C : \z\ < 1}- Since \\K\\ < 1, Lemma [2. 151 applies and yields that 
(/ - KK*) l l 2 K = K(I- K*K) l l 2 . Therefore, (EElQl) simplifies to nothing but the identity X = K. 
Now the claims follow from the inclusion (12.221) and the identities (12.191 ). (12-2 lb . □ 

3. Bounds on uniformly contractive solutions to the Riccati equations 

Assuming Assumption 12.11 in this section we prove several norm bounds on uniformly con- 
tractive solutions K to the Riccati equation (11.111 ) (provided such solutions exist). These bounds 
are obtained under the hypothesis that either the spectra of the operators Zq = Ao + BK and A i 
or the spectra of Zq and Z\ = A\— B*K* are disjoint. Note that, by Theorem 12. lOI (i). the as- 
sumption \\K\\ < 1 implies that the spectra of Zo and Z\ are both real, that is, spec(Zo) C R and 
spec(Z!) C R. 

Throughout this section we use the following notations: 

5z ,a, :=dist(spec(Z ),spec(A!)), (3.1) 
Szo,Zi :=dist(spec(Z ),spec(Zi)). (3.2) 

3.1. Semi-a posteriori bounds. First, we establish norm bounds on K that only contain the 
norm of B and the distance SzqAi- Therefore, these bounds may be viewed as semi-a posteriori 
estimates on K since the set spec(Zo) = spec(Ao) corresponds to the perturbed operator L = A + 
V (see Theorem 12 .10b . while the other set, spec(Ai), is part of the spectrum of the unperturbed 
operator A. 

Theorem 3.1. Assume that Assumption |Z71 holds and suppose that the Riccati equation (11.111 ) 
has a weak (and hence strong) solution K G ^(^1,^0) such that \\K\\ < 1. Then: 

(i) If the spectra of the operators Zq = Aq + BK, Dom(Zo) = Dom(Ao), and of A\ are 
disjoint, i.e. 

5zo Ai > 0, (3.3) 
then the solution K satisfies the inequality 

„ „ 71 IIBII 

\\K\\<--^-. (3.4) 

(ii) If, in addition, one of the sets spec(Zo) or spec (A 1) lies in a finite or infinite gap of the 
other one, i.e. 

conv(spec(Z )) nspec(Ai) = (3.5) 

or 

spec(Zo) nconv(spec(Ai)) =0, (3.6) 
then the solution K satisfies the stronger inequality 

„ „ l|5|| 

||*|| < (3.7) 
°Z Ai 



BOUNDS ON THE SPECTRUM AND REDUCING SUBSPACES 1 1 

Proof. The assumption that K is a strong solution to the Riccati equation (11.111) is equivalent to 

Ran(*| Dom(Zo) ) = Ran(*| Dom{Ao) ) C Dom(Ai) and 

KZqx-AiKx = -B*x for all x G Dom(A ) = Dom(Z ). (3.8) 

Since K is a uniform contraction, ||*|| < 1, we can use Theorem l2.10l (i) to rewrite (13.81 ) as 

K(I-K*Ky l/2 A (I- K*K) l/2 x-A l Kx = -B*x for all x G Dom(Ao), (3.9) 

where A is the self-adjoint operator defined by ( I2.9I ). If we choose x = (I — K*K)~ l l 2 y with 
y G Dom(Ao), we can write (13.91 ) as 

K(I - K*Ky l/2 A Q y-A l K(I - K*Ky l/2 y = -B* (I - K*Ky l/2 y (3.10) 

for all y G Dom(A ); 

note that Ran*(7 - ***r 1/2 | Dom(Ao) C Dom^j) since Ran(7 - ***) 1/2 | Dom(Ao) = Dom(A ) 
(see d2.9l >) and Ran(*| Dom ^) C Dom(Ai). Equality (13.101 ) means that the operator 

X = K(I -K*Ky l/2 (3.11) 
is a strong solution to the operator Sylvester equation 

XA Q - A Y X = Y (3.12) 
with Y = —B* (7 - K*K)~ X I 2 . Obviously, for the norm of Y we have 

II Dll 

llyll < H I' (3.13) 

If \K\ = \/K*K denotes the modulus of K, then the modulus |X| = y/X*X of the operator X 
denned in (13.111 ) is given by 

|X| = |*| 2 ) 1/2 . 

Taking into account that |||X||| = ||X|| and |||*||| = ||*||, the spectral theorem implies that 

||X|| = — J3=. (3.14) 

Due to the similarity ( 12.91 ) of the operators Ao and Zo, we have spec(Ao) = spec(Zo) and hence, 
by (fH, 

dist(spec(A ),spec(Ai)) =5^. (3.15) 
Applying Theorem 1 2 . 1 3 1 and using (13.131) as well as (13- 14b . we readily arrive at 

11*11 „ „ ||y|| iiflii 

X <c-£-±- <c- 



v i ~\m 2 5z °^ ~ 5 Zo4i \/i-ii*ii 2 

where c = n/2 in case (i) and c = 1 in case (ii) so that, in both cases, 



\K\\<c-f^-. □ 



Remark 3.2. In order to compete with the hypothesis ||*|| < 1, the bounds (13.41 ) and (13.71 ) are 
of interest only if ||S|| < 25z oAl /7T in case (i) and ||S|| < 5z ,Ai in case (ii)- 

Remark 3.3. For all spectral dispositions such that (13.51 ) or (13.61) holds and ||S|| < 5z 0jAl , the 
bound (13.71) is sharp in the sense that given an arbitrary /3 > and arbitrary 8 > j3 one can 
always find A and V such that ||V|| = /J, 5z .a, = 8, and ||*|| = fi/S. 



12 S. ALBEVERIO, A. K. MOTOVILOV, AND C. TRETTER 

The following examples illustrate the sharpness of (13.71) and Remark [331 
Example 3.4. Let Sjo = C 2 and Sji = C. Assume that b, d € M are such that < b < d/2 and let 

d J' Al=0 ' B ={b 

For this choice of Aq, A i , and B, the Riccati equation (11.111 ) has a unique uniformly contractive 
solution of the form K = (0 — Kj where K is given by 

h 

(3.16) 



Hence 



o f + Vf-^ 2 



and the set spec(Ai) = {0} lies in the gap (—d,d/2 + ^/d 2 /A — b 2 ) of spec(Zo), so that (13.61 ) 
holds. Altogether we have 

„ „ l|5|| 

||^|| = iUL, (3.17) 

i.e. equality in (13.7I ). 

Example 3.5. Let $j = C and Sji = C 2 . Assume that b,d £ M are such that < b < d and set 

- d °\ S=( fo fo 
By inspection, one can verify that the 2x1 matrix 



b 

\fld 
b 

s/2d 



K = 

is a solution to the operator Riccati equation (11-1 lb - Clearly, 

n „ „ „ b 

\\B\\=b, \\K\\ = - (3.18) 



and 

Z =A + BK = 0, Z l =A l -B*K*=[ a ^ 2* . 

\ Id " 2d ) 

Obviously, the set spec(Zo) = {0} lies within the gap (—d,d) of the set spec(Ai) = {— d,d}. 
Furthermore, Sz^, = d and hence, by (13.18b . 

||*|| = -jp!L (3.19) 
°z Ai 



For later reference, we note that spec(Zi) = { — Vd 2 — b 2 , Vd 2 — b 2 } so that 5z 0i Zi = Vd 2 — b 2 
and thus 

, „ \\B\\ \\B\\ 

\k\\ = , y y = y 11 0.20) 

Sl Zl +b 2 J8l Zi + \\B\\ 2 
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Example 3.6. Let fjo = fii = C. Assume that b,d G R are such that < b < d/2 and set 

A Q = -d/2, Ax = d/2, B = b. 

Then the Riccati equation (11.111 ) appears to be the numeric quadratic equation bK 2 + Kd = — b. 
The only solution * = k G R with norm ||*|| = < 1 where fc is again given by (13.16b . One 
immediately verifies that 



Z =A + BK = -^-b 2 , Z X =A,-B*K* = ] Jf--b 2 . (3.21) 

Here the sets spec(Zo) and spec (A i) are even subordinated to each other, so that both (13.51 ) and 
(13.61 ) hold. Together with ( 13.161 ) and ( 13.211 ), we again obtain the equality 

||tf|| = -pi. (3.22) 

°Z ,A l 



For later reference, we also observe that 

11*11 \\B\ 



l-||*ll 2 5; 



(3.23) 



Zo,Zi 



3.2. Completely a posteriori bounds. In this subsection we consider the case where the spec- 
tra of the operators Zo = Aq + BK and Z\ =A\— B*K* are disjoint. The bounds on * obtained 
here depend only on ||B|| and on the distance 8z Q ,z l between the subsets spec(Z ) = spec(A ) and 
spec(Zi) = spec(Ai) of the spectrum of the perturbed operator L = A + V (see Theorem 12. 101) . 
Therefore, they may be viewed as a posteriori bounds on * 

We begin with the most general result where nothing is known on the mutual position of 
spec(Zo) and spec(Zi) except that they do not intersect. 

Theorem 3.7. Assume that L = A + V satisfies Assumption \2.1\ and suppose that the Riccati 
equation (11.111 ) has a weak (and hence strong) solution * € 3$(S}i,$)o) such that \\K\\ < 1. 
If the spectra of the operators Zq = Ao + BK, Dom(Zo) = Dom(Ai), and Z\ = A\ — B*K*, 
Dom(Zi) = Dom(Ai), do not intersect, that is, 

Szo,Zi > 0, 

then 

\\K\\ < --JUL. (3.24) 

Proof. By Lemma 12.161 the Riccati equation (11.111 ) can be written in the form (12.161) . For the 
term Y = -(I - KK*) 1 I 2 B*{1 - K*K) l l 2 on the right-hand side of (ETTBl we have ||Y|| < ||B|| 
since both *** and *** are non-negative and, in addition, ||*||||**|| < 1. Since spec(Ao) = 
spec(Zo) andspec(Ai) = spec(Zi), we have dist(spec(Ao),spec(Ai)) =5z ,z 1 - To complete the 
proof, it remains to apply the bound (12.111) from Theorem l2.13l to (12.161) . □ 

Remark 3.8. Under the stronger assumption that one of the spectral sets spec(Zo) and spec(Zi) 
lies in a finite or infinite gap of the other one, i.e. if 

conv(spec(Z,)) n spec(Zi_ ; ) = for some i = 0, 1, 

Theorem 12. 1 3 1 also yields the estimate 

ll*ll < Jul. 0.25) 

d Zo,Z! 
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This estimate, however, is of no interest in the case where the corresponding operator Z, is 
bounded: the bound (13.271) in the following theorem is stronger than (13.25I ). 

Theorem 3.9. Assume that L = A+V satisfies Assumption \2.1\ and suppose, in addition, that 
Aq is bounded. Let the Riccati equation (11.111 ) have a weak (and hence strong) solution K £ 
3&(f)i,Sjo) such that \\K\\ < 1 . If the spectrum of the ( bounded) operator Zq = Aq + BK lies in 
a finite or infinite gap of the spectrum of the operator Z\ =A\ —B*K*, Dom(Zi) = Dom(Ai), 
that is, 

conv(spec(Zo)) n spec(Zi) = 0, (3.26) 

then 

\\k\\ < y y =. (3.27) 
+ II*" 2 

Proof. Throughout the proof we assume that B ^ and, thus, necessarily 

K^O. (3.28) 

Let U be the partial isometry in the polar decomposition K = U\K\ of K. If we adopt the 
convention that U is extended to Ker(X) = Ker(|^T|) by 

U\K S r(K)=0, (3.29) 

then U is uniquely defined on the whole space f)o (see [ 10, Theorem 8. 1.2] or |fT7l §VI.7.2]) and 

U is an isometry on Ran(|^|) = Ran(^*). (3.30) 

First we apply Lemma 12.161 and transform the Riccati equation (11.11b to the form (12.161) . 
Since the operator Ao is bounded, Ao G 38{ff)o), we may then rewrite (12.161 ) as 

KAq-MU\K\ = -(I-KK*) l/2 B*(I-\K\ 2 ) l/2 = -B*(I- \K\ 2 ) 1/2 (3.31) 

where we have set 

B = B(I-KK*) 1/2 . (3.32) 
We tackle the cases where the gap of spec(Ai ) containing the set spec(Ao) is finite or infinite 
in a slightly different manner. If this gap is finite, we may assume without loss of generality that 
it is centered at zero, i.e. it is of the form (—a, a) with a > 0; otherwise, we simply replace Ao 
and A] in (13.311 ) by A = Ao — X\I and Aj = Ai — X\I, respectively, where X\ is the center of 
the gap. Then 

Oep(Ai), IIA^ 1 !! < -, and ||A || <a-8 Zo ,z r (3.33) 

If the gap of spec(Ai) containing the spec(A ) is infinite, we may assume without loss of gener- 
ality that the interval [min spec (Ao), max spec (Ao)] is centered at zero and that the spectrum of 
Ai lies either in the interval (— °o ; — a] where a = — max spec (Ai) or in the interval [a,°°) where 
a = min spec (Ai). Then, again, all three statements of (13.331 ) hold. 

In the following, we may thus treat the two above cases together. Since spec(Ai), we 
further rewrite (13.311) in the form 

U\K\ =A~ l (kA + B*(I-\K\ 2 ) 1/2 ) (3.34) 



and set K = max spec(|A"|) = \\K\\. By assumption (13.281 ). we have fc> 0. 

If fc is an eigenvalue of \K\ and x a corresponding eigenvector with ||x|| = 1, then, by applying 
both sides of (13.341 ) to x, we immediately arrive at 



kUx = A' 1 (KAox + \Jl - K 2 B*x) . (3.35) 
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By (El) , we have < < ||5||||(/-Z'Jr , ') 1 / 2 || < Then, by ([333]) and 035]), we 

obtain 

jc< i(jc(a-5 ZojZl ) + Vl- K- 2 ||S||), (3.36) 

taking into account that x G Ran(|7v|) and thus \\Ux\\ = \\x\\ = 1 by (13.301 ). 

If K is not an eigenvalue of \K\, then it belongs to the essential spectrum of \K\. Hence we 
obtain a singular sequence {x n }™ =i of \K\ at K"by choosing arbitrary 



x„eRanE|jq((ir(l-l/n),ir]), \\x„\\ = 1, n = l,2,...; 



(3.37) 



here E\ K \ denotes the spectral measure of |^T| and, at the same time, the (right-continuous) spec- 
tral function of |7sT|, that is, E|jf|(ju) = E^ ((— °°,ju]). Obviously, \K\x n = Kx n + £„ with 



Similarly, (7 - \K\ 2 ) l / 2 x n = y/l - K 2 x n + £, with 



dE\ K \ (/x)(/x-K)x n 



< 



(3.38) 



Knl 



< sup 

/xe(K(i-i/»),K] 



1 - fC Z 1 



v/l-/i 2 -\/l-K 2 | 



2\ ! /2 



(1 - rc 2 ) 1 / 2 



< 



Applying both sides of equality (13.341) to x n , we arrive at 



KUx n = A, 1 (KA x n + \/l-K 2 B*x n ) + j3„ 



with 



j3„ = A^TT - £/£„ as n -> °°; 



here we have used that e„ = |*T|x„ — fCtn — ► and = (7 — ITy) 2 ) 1 / 2 *;,, — V 1 — (Ac,, 
n — > oo by (13.381 ) and ( 13.391 ), respectively. Since x„ G Ran(|7v|) and thus ||f/;c n || = ||x„|| 
(13.301) . the relation (13.401 ) implies that 

K<- {K(a- Sz, A ) + Vl ~ K 2 \\B\\) + Hft.ll, «=l,2,...,oo, 



(3.39) 



(3.40) 



(3.41) 

-> as 
= 1 by 



(3.42) 



which, by (13-4 lb . turns into the bound (13.361 ) after taking the limit n — > oo. Solving inequal- 
ity (13.361) for K and recalling that K = \\K\\,v/e conclude the estimate (13.27I ). □ 

Remark 3.10. The bound (13.271) is sharp. In fact, equality (13.201) in Example |3 .5 1 shows that, for 
the spectral dispositions (13.261 ) where spec(Zo) lies in a finite gap of spec(Zi), equality prevails 
in (13371 . 

The strongest a posteriori bound for the solution K is obtained under the assumptions that the 
spectra of Zq and Z\ are subordinated, i.e. 



max spec(Zo) < min spec(Zi) or max spec(Zi) < min spec(Zo), 
and that Aq and A i are bounded. 



(3.43) 



16 



S. ALBEVERIO, A. K. MOTOVILOV, AND C. TRETTER 



Theorem 3.11. Assume that L = A + V satisfies Assumption 12. 1 1 and suppose, in addition, that 
the operators Aq, A\ (and hence A and L) are bounded. Let the Riccati equation (11.111 ) have 
a weak (and hence strong) solution K G &($)i,Sjo) such that \\K\\ < 1. If the spectra of the 
operators Zq = Aq + BK and Z\ = A \ — B*K* are subordinated, that is, 

max spec (Zo) < min spec (Zi) or max spec (Zi) < min spec (Zq), (3.44) 

then 

\\K\\ < tan - arctan — ^ — ^- . (3.45) 
V 2 °z .zj 

Proof. As in the proof of Theorem 13.91 we apply Lemma |2.16| and rewrite the Riccati equation 
in the form (12.161) ; note that the self-adjoint operators Ao and Ai on the left-hand side of (12.161 ) 
are bounded by Theorem 12.101 since Ao and A i are bounded, B is bounded, and K is a uniform 
contraction. Assuming that B ^ 0, we again have K / (cf. (I3.28I )). 

Let U be the partial isometry in the polar decomposition K = U\K\ of K (see the proof of 
Theorem[L9]>- By Lemma [245] with <p(z) = (1 -z) 1/2 , M = U*, and N = \K\ 2 U*, we obtain 

U*(I- KK* ) 1/2 = U*(I-U \K\ 2 U* ) 1/2 

= (I-U*U \K\ 2 ) l/2 U* 

= (I-\K\ 2 y/ 2 U*. 

Here, in the last step, we have used the property that U is an isometry on Ran(|i<T|) = Ran(/r*) 
by (13301 so that U*U\K\ = \K\ and U*U\K\ 2 = \K\ 2 . 

If we apply the operator U* to both sides of (12.161 ) from the left, we arrive at an equation that 
only involves \K\, but not K and K* themselves: 

\K\Ao-U*AiU\K\ = -(/- |#| 2 ) 1/2 £/*fi*(/-|^| 2 ) 1/2 . (3.46) 

Let K = max spec(|^|). Clearly, < K = \\K\\ < 1. If K" is an eigenvalue of and x G Sjo, 
\\x\\ = 1, is an eigenvector of at K, that is, \K\x = Kx, then (13.461 ) immediately implies that 

— ^((A x,x) - (AxUx^Ux)) = -(U*B*x,x). (3.47) 

Since x G Ran(|^|), by ( 13.301 ) we have \\Ux\\ = \\x\\ = 1 so that 

(AiUx,Ux) > min spec(Ai) = min spec (Z\ ) , (3.48) 

(Aqx,x) < max spec(Ao) = max spec(Zo). (3.49) 

Because the spectra of Zq and Z\ are subordinated by assumption ( 13.441 ), the inequalities d3.48| ) 
and ( 13.491 ) yield that 

|(A x,x)-(A 1 [/x,^)| >5 Zo , Zl >0. 
This and ( 13.471 ) imply the inequality 

K \\B\\ 

—— 2 < J-JL. (3.50) 

1 - K dzo.Z, 

If K is not an eigenvalue of \K\, it belongs to the essential spectrum of \K\. We introduce a 
singular sequence {x„}~ =1 of \K\ at K as in (13.371) ; in particular, ||x„|| = 1. With this choice 
of x„, (13.461 ) implies that 

1 {(AQX n ,x n )-(A x Ux m Ux n )) = -(U*B*x n ,x n ) + a n (3.51) 
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where 

_ (U*A l Ue n ,x n )-(A x n ,e n )-(U*B*i; n ,i; n ) (U*B*£ n ,x n ) + (U*B*x n ,£ tt ) 

Because of x n G Ran(|X"|) and (I3.30I ). we have ||£/x„|| = \\x n \\ = 1. Now the same reasoning as 
in (137481 ) and tf3~49T> yields that 

(A x„,x n )-(A l Ux n ,Ux n )\ > 5 Zo , Zl >0. 

Hence (13.511) shows that 



LB cc n 
i_n_ + 

3Z 0! Z, 5 Zo! z, 



■7— 2<^ + ^, «=l,2,...,oo. (3.52) 

1 TC 



As ||jc„ || = 1 and both Ao and Ai are bounded operators, (13.381 ) and ( 13.391 ) show that a n — > for 
n — > oo. Taking the limit « — > oo in (13.521 ). we again arrive at inequality (13.501) . 

To complete the proof it remains to notice that, by the formula for double arguments of the 
tangent function, the left-hand side of (13.501 ) may be written as 

K 1 

j = ~ tan(2arctan fc), (3.53) 

and to recall that fc = \\K\\. □ 

Remark 3.12. The bound (13.451 ) is optimal. This may be seen from Example 13.61 where Sjo = 
Sji =CandA = -d/2,A\ = d/2, B = b with b,d G K, < b < d/2. In fact, by 0331 , equality 
( 13.231 ) therein is equivalent to 

,1 2\\B 
\K\\ = tan - arctan 



2 5 Zo,Z! 

4. Proofs of Theorems II .21 and II .31 

Using the results of Section |3l we are now able to prove our main results, Theorems 1 1.21 and 
Theorem 1 1.31 which were formulated in the introduction. In particular, Theorem 11.21 appears to 
be a corollary to Theorem l3.1l in Section[3] 

Proof of Theorem lL2\ According to the definitions (11.8I ). (11.9b and Theorem l2.10l we have 

Oi = spec (Aj), ol = spec (Z;) , i = 0, 1 , 

with Z = A + BK, Dom(Z ) = Dom(A ) and Z\ =A\ — B*K* , Dom(Zi ) = Dom(A i ) as above. 

We prove the theorem in the case dist(a ( ), G\) = SzqAi > 0; the case dist(a 1 / , ao) > may be 
reduced to the first case by replacing the involution / with /' = — J and making the corresponding 
index changes in the notations of Assumption 12. II 

By assumption, fj is a maximal uniformly positive subspace of the Krein space 8. = {ft, J}. 
Thus Lemma [2781 implies that f) is the graph of a uniform contraction K : f)o — ► -fh, i-e. f) = 
&(K). By assumption, S)' is also a reducing and hence invariant subspace of L. Now Lemma 
I2.7l yields that K is a strong solution to the operator Riccati equation (11.1 lb . By Lemma 12791 and 
Theorem EjOl we know that Sj\ = ( S(K*)= &(K)^ = and hence, by 423} and definition 

dunii, 

||tan0 o || =0(£o,#o) = ll^ll = \\Z*\\=®(fiiM = ||tan©i||. (4.1) 
Now both claims (i) and (ii) are immediate consequences of the respective statements (i) and (ii) 
of Theorem 13. II □ 
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Remark 4.1. Under the natural assumption that || V \\ < 5, for some i G {0, 1}, the bound dl.191 ) is 
optimal with respect to the mutual positions of the sets a, and o[_ i described in condition (1 1 - 18b - 
This follows from Remark [331 and the subsequent examples together with the equalities d4.1| ). 

Theorems 13.71 13.91 and [3_7TT] enable us to prove Theorem [T3] 

Proof of Theorem 1731 By definition dL9b and Theorem l2.101 we have 

a/ = spec(Z;), i = 0,l, 

with Zo = Ao + BK, Dom(Zo) = Dom(A ) and Zi = Ai - B*K*, Dom(Zi) =Dom(Ai). Hence 
assumption (11.201 ) implies that 5z .Zi = dist(a ( ), o[) = 8 > by (13.2I ). 

As in the proof of Theorem 11.21 we conclude that fj is the graph ^(K) of a uniformly 
contractive strong solution K to the operator Riccati equation (ll.lll ). while 9)\ is the graph 
&(K*) of the adjoint of K. 

For claim (i), the bound (11.211 ) follows from estimate (13.241 ) in Theorem 13.71 using relation 
CD. 

For claim (ii), the bound (11.23b for i = follows from estimate (13.271 ) in Theorem |3.91 again 
using relation (14.1b ; for i = 1 it follows from the case i = by passing from / to the new 
involution /' = — /. 

For claim (iii), the bound (11.25b follows from estimate (13.45b in Theorem l3.11l if we use (14.1b 
and the facts that ||tan@ y || =tan||0 y || and, by (11.12b . ||tan20 ; -|| = tan2||@y||, j = 0, 1. □ 

Remark 4.2. The bound (11.23b is sharp whenever the gap of a' l _ j containing o[ is finite. For 
i = this follows from the fact that the estimate (13.27b in Theorem 13.91 is sharp by Example 
13.51 (see Remark [3 .10b together with the identity (14.1b : for i = 1 it follows using the involution 
f = —J instead of J. 

Remark 4.3. The bound (11.25b is best possible. This follows from the fact that the estimate 
(13.45b in Theorem 13.1 II is sharp by Example 13.61 (see Remark 13.12b together with (11.12b and 

CD. 

5. Estimates for the perturbed spectra 

In the next section we want to use the operator angle bounds of Theorems 11.21 and 11.31 to 
prove a priori bounds on the variation of spectral subspaces of the self-adjoint operator A under 
an off-diagonal /-self-adjoint perturbation V. To this end, we establish some tight enclosures 
for the spectral components of the perturbed operator L = A + V in the present section. 

We assume that the initial spectra Go = spec(Ao) and Oy = spec (A i) of the block diagonal 
entries Aq and A\ of A (see (11.4b ) are disjoint, i.e. 

d = dist( Ob, di) >0. (5.1) 

Then the subspaces f)o and f) i introduced in Assumption 12.11 are the spectral subspaces of A 
associated with the spectral components Oq and 0\ , respectively. 

In the following, our aim is to find certain bounds for the perturbation V and a constant ry > 
such that 

dist(c^, o[ ) > and a/ C O rv (a,) nM, £ = 0,1. (5.2) 
This yields the lower bounds 

8i = dist(a,-,cj_,-) > d-r v , i = 0,l, and 8 = dist(ag,a 1 / ) > d-2r v . (5.3) 
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Together with the estimates in Theorems 1 1.21 and 1 1.3 1 respectively, they will give us the desired 
a priori estimates for tan© ; , depending only on the initial distance d of the unperturbed spectra 
and on the norm of V. 

For completely arbitrary (i.e. not necessarily off-diagonal) perturbations V of the self-adjoint 
operator A, it is well-known that the assumption 

IM| < d - (5.4) 

guarantees that (15.21 ) holds with 

ry = ||V|| (5.5) 
(see, e.g., [I7j Section V.4]) and hence, by (15.31 ). 

5i>d-\\V\\ > |, i = 0,l, 8 >d-2||V|| >0. (5.6) 

For off-diagonal perturbations V, earlier results in |[T9l . ll20l for self-adjoint V and in |6l for 
non-symmetric V show that the constant ry may be improved considerably. In the following we 
extend these results under the sole assumption d5.ll ) that the spectral components ct an d G\ of 
A are disjoint. 

Unlike the previous sections, we do not always require A to be self-adjoint and V to be /-self- 
adjoint; here we use the following more general setting. 

Assumption 5.1. Let Sjq and S)i be complementary orthogonal subspaces of the Hilbert space Sj. 
Assume that A is a closed operator on fj = S)o © S)i diagonal with respect to this decomposi- 
tion, i.e. 

' An 



^° jM, Dom(A) =Dom(A )©Dom(Ai), 



where Aq and A \ are closed operators on S)q and S)q, respectively. Suppose that V is an off- 
diagonal bounded operator on Sj, i.e. 

v=(° B 

\ C 

withBe &{$)i,fio), C G @($)o,S)i), andlet 

L = A + V=(^ A ® B Ax \ Dom(L) = Dom(A). (5.7) 

The following two elementary auxiliary results are used in the proofs below. 

Lemma 5.2. Assume that L=A + V satisfies Assumption \5. 1 1 and define the Schur complement 
So of A by 

So (A) :=Aq-X-B(Ai -iy l C, Dom(S (A)) :=Dom(A ), 
for A G p(Ai). Then 

(i) the resolvent set p {Sq) := {A G p(Ai) : So (A) is bijective] of So satisfies 

p(So)=p(L)np(A l ); 

(ii) for A G p (A ) D p (Ai ) we have 

\\B(Ai -A)" 1 C(A -A)" 1 || < 1 => Aep(L). (5.8) 
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Proof. Both claims (i) and (ii) are well known (see, e.g., |[32l ): we recall the short proofs for the 
convenience of the reader. 

(i) It is easy to check that, for arbitrary / G S)o, g G $)i and x G Dom(Ao), y G Dom(Ai), 



(L-A) 



yyj \g) ^ \ y =(A 1 -X)-\g-Cx), 

which proves the claim. 

(ii) For A G p(Aq) Hp(A\), one can write 

So (A) = {l—B(A\ - X)- l C(A - A)" 1 ) (A - A). 

Thus a Neumann series argument together with (i) proves (15 - 8b - □ 

Lemma 5.3. Lef a, £>, v G M. be such that 8 := b — a > ant/ < v < 8/2. Then 

(t -a)(b-t) > v 2 <==^ a + r<t <b-r 

where 



8 S 2 i /I . 2v\ 

r=-- W — -v 2 = vtanl -arcsin— I . (5.9) 

Proof. The claims are obvious; for the last equality, observe the formula for double arguments 
of the sine function in terms of the tangent function. □ 

In the following theorem, we consider the case where the diagonal entries Aq and A\ of 
the block operator matrix (15.71) are self-adjoint and their spectra do not intersect; the bounded 
perturbation V need not have any symmetry here. 

Theorem 5.4. Assume that L = A + V satisfies Assumption ^, h and let Aq andA \ be self-adjoint. 
Assume, in addition, that their spectra (To = spec(Ao) and di = spec (A i) are disjoint, i.e. 

d = dist(cb,ai) > 0, 

and let the entries B and CofVbe such that 



Then 

spec(L) = ad 0a[, a[ C 0,- v (a,-) n R, i = 0, 1, (5. 1 1) 

where ry is given by 



r v 



y/\\B\\ \\C\\ tan ^arcsin ^^ H j < v/||5|| ||C||; (5.12) 

in particular, if V is J -self -adjoint and hence C = —B*, then 

/ 1 2||V||\ 

r v = \\V\\ tan ( -arcsin-^- 11 j < ||V||. (5.13) 

Proof. Throughout the proof, we assume that A G C is such that 

dist(A,a Uai) > r v ; (5.14) 
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hence, in particular, A G p(Ao) Hp(Ai). Since Aq and Ai are assumed to be self-adjoint, we 
have ||(A ; - = l/dist(A, a,-), i = 0,l, and thus 

[^(Ai-Ar^Ao-A)- 1 !!^ - M^jL -r . (5-15) 

dist(A,ao)dist(A,ai) 

First we consider the case that A lies in a strip of the form 

{zeC\a + r v < Rez <b-r v } (5.16) 

where (a, ft) is a finite gap of the spectrum of A with a G Go and ft G <n; the case a G di and 
ft G Ob is analogous. Then we have b — a>d and hence, by assumption ( 15.101 ) and Lemma [531 
we obtain 

llflllllCH IIBIIIICH l|fi||||C|| 

< .„ . " 1 ' " . . < 1. (5.17) 



dist(A,CT )dist(A,CTi) \X-a\\X-b\ ~ (Re A — a)(ft — ReA) 

Now (15.151 ) and (15.171 ) together with Lemma E2](ii) show that A G p(L). 

If A does not belong to a strip of the form (15.161 ) with a G a, and ft G (7i_; for z = or i = 1, 
then it is not difficult to check that, by (15 - 14b . either for i = or i = 1 

dist(ReA,a,) > d-r v , dist(A,<Ti_,) > ry. (5.18) 

Combining (15- 15b . (15.181 ) together with the definition (15.121 ) of ry, we arrive at the estimate 

\\B(Ai - A)- 1 C(A -A)- 1 || < l|g ] ll|C|1 . = 1. 

r v {d-r v ) 

Hence Lemma l5T2l (ii) again shows that A G p(L). □ 

Remark 5.5. Theorem [53] improves the spectral bounds given in (6l Remarks 4.6 and 4.13]. In 
fact, the bound r in J6l (4.1 1)] can be written equivalently as 

r = ^/pfljcf tanh ^ arctanh ^ZHH^ = ^/||B|| ||c|| tan ^ arcsin ^ZHH^ (5. 19) 

and applies whenever ||fi||||C|| < 8/2. Here, by O Theorem 4.11], we have 8 = 2d/n < d in 
the general case (15.11 ) and 5 = d if one additionally assumes that 

conv(a )nai =0 or a nconv(ai) = 0. (5.20) 



This shows that, under the additional assumption (15.201) . the spectral bound in Theorem 
coincides with the one in [6, Remarks 4.6], whereas in the general case d5.ll ) Theorem |5.4| holds 
with the weaker norm bound ||C|| < d/2 on V and the bound ry is strictly smaller than the 
bound r in |6j (4.11)]. 

Remark 5.6. The spectral bound (15.1 II ) with ry given by (15.121 ) is optimal. This may be seen 
from Examples 4.15, 4.16]. 

In the next two theorems we drop the assumption that A is self-adjoint. Instead we impose 
conditions to ensure that the components Ao and Ai of A satisfy certain resolvent estimates. 

To this end, we use the numerical range W(T) of a linear operator T with domain Dom(r) 
in a Hilbert space, defined as 

W(T) = {(Tx,x) I x G Dom(T), ||jc|| = 1}. 
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Recall that the numerical range is always convex and that spec(r) C W(T) if every (of the at 
most two) connected components of C\W(T) contains at least one point of p(T) (see |[T71 
Theorems V.3.1 and V.3.2]); in this case, 



'" r - A »'''^ d. s ,(A,W» ' 
First we consider the case where the spectra and the numerical ranges of Aq and A\ are 
separated by a vertical strip. 

Theorem 5.7. Assume that L = A + V satisfies Assumption 15. H and that there exist a,b Gl such 
that d = b — a > 0, 

ReW(Ao) <a<Z><ReW(Ai) or Re W(Ai) < a < b < Re W(Aq), (5.21) 

and 

{z.£C\a< Rez. < b) n p(A ) n p(Ai) / 0. (5.22) 



VMM < d/2, 

and ry is defined as in ( 15.121 ), then 

{z G C | a + r v < Rez <b-r v } C p(L). 

Proof. Without loss of generality, we suppose that ReW(Ao) <a<b< ReW(Ai). In this case 
the assumptions (15.211) . (15.221 ) imply that (see 1 17. Theorem V.3.2]) 



spec(A ) C W(A ) C {z € C | Rez < a}, spec(Ai) c W(A r ) C {z € C | Rez > b}. (5.23) 
LetA GCbe suchthatc? + r y <ReA <b + r v and hence A G p(A )np(Ai). Then, by (15.231 ). 

Il^o-A)- 1 !^— 7 -i-_<— i_ , ||(Ai - A)' 1 )) < , 1 < 



dist(A,W(A )) ~ ReA-a' " v ; " ~ dist (X,W(Ai)) ~ b-ReX ' 
Thus Lemma [531 shows that 

|| B (A,-^C(A -^||< J1M_ <1 , 

and hence A G p (L) by Lemma l5T2l (ii). □ 

Next we consider the case where the spectra and the numerical ranges of Ao and A i (and 
hence of A) lie in one half -plane and the perturbation V is /-self-adjoint, i.e. C = — B* . 

While all previous theorems were of perturbational character, the following theorem is not. 
In fact, we prove implications of the form 

Re spec(A) < ReW(A) < a =>• Re spec(L) < a, 

independently of the norm of V. 

This type of results relies on the quadratic numerical range W 2 (L) of the operator L with 
respect to the block representation (15.7I ). The set W 2 (L) is defined as (see l25l (2.2)] and also 
Il33l Definition 3.1] or OH Definition 2.5.1]) 

W 2 (L) = |J spec(L, !V ) 

A'GDom(A()). vGDom(A^ ) 

IWHMhi 
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where L x>y G M2(C) is a 2 x 2 matrix given by 

^ '= ( %°x,y) (My*y) ) ' X ^ Dom ( A °)^ G Dom(A 1 ), ||*|| = \\y\\ = 1. 

The quadratic numerical range is not convex and may consist of at most two connected 
components. It is always contained in the numerical range, W 2 (L) CW(L), and the inclusion 
spec(L) C W 2 (L) holds if every connected component of C\ W 2 (L) contains at least one point of 
p(L) (see (25] Theorem 2.1], (33J Proposition 3.2, Theorems 4.2 and 4.7] or (32j Theorem 2.5.3, 
Theorems 2.5.10 and 2.5.15])). 

Theorem 5.8. Assume that L = A + V satisfies Assumption 15. 1 \ with C = —B* and let a,b 6l 

(i) // ReW(A) < a and p(A) n {z G 

(ii) // Re W (A) > b and p(A) n{z£ 
in particular, if A is self-adjoint, we have 



Rez > a] 7^ 0, then Re spec(L) < a. 
Rez < b) 7^ 0, f/jerc Re spec(L) > Z>; 



inf spec (A) < Re spec(L) < sup spec (A). 

Proof. We prove (i); the proof of (ii) is completely analogous. By the assumption on A and 
since V is bounded, it is obvious that |zGC| Re z > np(L) ^ for L = A + V. Hence 
spec(L) C W (L) and so it suffices to show that ReW 2 (L) < a. 

Suppose, to the contrary, that there exists a A £ W 2 (L) with ReA > a. By the definition of 
W 2 (L), there are x G Dom(Ao), y G Dom(Ai), ||jc|| = \\y\\ = 1, such that X G spec(L vy ), i.e. 

= det (L r , y - X) = ((Aox,x) - A) ((A iy ,y) - A) + |(Bj,x)| 2 , 

where we have used that C = — B*. Splitting into real and imaginary parts, we conclude that 

= Re ((A Q x,x) - A) Im ((A iy ,y) - A) +Im ((Aox,x) - A) Re ((A iy ,y) - A) , 

-\(By,x)\ 2 = Re {(A Q x,x) - A) Re ((A iy ,y) - A) - Im ((A x,x) - A) Im {(A iy ,y) - A) . 

Solving the first equation for Im ((A\y,y) — A) and inserting into the second equation, we find 

- \(By,x)\ 2 = ((Re(A x,x) -ReA) 2 + (lm(A x,x) -ImA) 2 ) ^^^Z^ , (5.24) 

By the assumption on A, we have Re A > a > ReW(A) = Re (W(A )UW(Ai)) and hence both 
the first and the second factor on the right hand side of (15.241 ) are positive, a contradiction. □ 



6. A PRIORI BOUNDS ON VARIATION OF SPECTRAL SUBSPACES 

In this section we use the (semi-) a posteriori norm bounds for the operator angles from 
Theorems 11.21 and 11.31 together with the spectral estimates from Section [5] to derive a priori 
estimates for the variation of the spectral subspaces of the self-adjoint operator A under a /-self- 
adjoint off-diagonal perturbation V. 

To ensure that solutions of the corresponding Riccati equations exist, we use some results of 
j6ll and ll35l . They provide sufficient conditions on the perturbation V and spectral sets Go and 
(7i guaranteeing that the perturbed operator L = A + V is similar to a self-adjoint operator on a 
Hilbert space and that the spectral subspaces of L associated with the perturbed spectral sets Oq 
and a[, both being real in this case, are maximal uniformly definite in the Krein space 
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In order to formulate the conditions from 051 . we need to specify all those finite gaps of the 
spectrum of A that separate the subsets Go and Of. We denote these gaps by A„, where n£Z 
runs from AL through N + with — oo < < 0, < < +oo, and let 

A„ = (a n ,b n ), -oo < ... < c?„_i < b n - X < a n < b n < a n+ \ < b n+ \ <...<°°, (6.1) 

assigning the value of n = 0, say, to the gap that is closest to the origin z = 0. For every such 
gap, we have A„ n Oq = 0, A„ n Of = 0, and a n G a,-, b n G a"i_j where either / = or / = 1. 
If the total number N of the gaps between Oq and 0\ is finite, then both AL, N + are finite and 
N = \N- 1 + N + + 1. Otherwise, at least one of N- and N + is infinite. 

The next theorem is an immediate consequence of (6J Theorem 5.8] and IT351 Theorem 3 and 
Corollary 4], combined with Theorem 15 .41 

Theorem 6.1. Assume that L=A + V satisfies Assumption \2.1\ and let the spectra Oq = spec (Aq) 
and <J[ = spec (A i) be disjoint, i.e. 

d = dist(cb,ai) > 0. 
Assume, in addition, that one of the following holds: 




Then 

spec (L) C Oq U o[ , o\ CO rv (Oi)f]R, i = 0, 1 , 

with ry given by ( I5.19I ). The operator L is similar to a self-adjoint operator on Sj. The spectral 
subspaces S)' and S)\ of L associated with the sets o' and o[ are mutually orthogonal in the 
Krein space {ft, J} and maximal uniformly positive resp. negative therein. 

Remark 6.2. In case (i) Theorem 16.11 follows from (6l Theorem 5.8 (i) and Remark 5.11], to- 
gether with Theorem 15.41 In case (ii) the spectral projections El(g-) of L associated with its 
isolated spectral components a/, i = 0, 1, are well-defined (see |[35l Corollary 4]), the operator 
J(El(Cq) — El(o{)) is self-adjoint, and there exists a ye (0, 1] such that 

7(E L ((j^-E L ((j()) >yl. (6.2) 

The latter was established (under much more general assumptions on V than (15.41) ) in the proof 
of |[35l Theorems 1 and 3]. Using inequality (16.21 ). one easily verifies that the spectral subspaces 
Sj'q = RanE^(aQ) and Sj\ = RanE^dj') are maximal uniformly positive and maximal uniformly 
negative, respectively. In fact, it suffices to show the uniform definiteness of one spectral sub- 
space (see Corollary 12.121) . 

Remark 6.3. The lengths of the gaps (a n ,b n ) of spec (A) separating the sets Go and Oi have to 
be uniformly bounded from below. Apart from this, condition (i) imposes no further restriction 
on the behaviour of the lengths, whereas condition (ii) requires that b n — a n tends to °° faster 
than \n\ as \n\ — > oo. 



The following a priori bound on the operator angles 0(f),-,i]-) between the unperturbed and 
the perturbed spectral subspaces Sjj of A and S) • of L = A + V improves the corresponding bound 
derived in [6, Theorem 5.8 (i)] (see Remark loTol below). 
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Theorem 6.4. Suppose that L = A + V satisfies the assumptions of Theorem I6.il Let £),-, 5j • 
be the spectral subspaces of A corresponding to a, and of L corresponding to o\, respectively, 
i = 0, 1. Then the operator angles &j = 0(f),-, f)-), i = 0, 1, satisfy the estimate 



K (\ 2||V| 

tan©, < — tan - arcsin ) , i = 0, 1: (6.3) 

2 V2 d 



conv(a )nai = or a nconv(ai) = 0, (6.4) 
/I 2||V||\ 

tan©,- < tan - arcsin , i = 0, 1. (6.5) 



$ in addition, 
then 



Proof. By Theorem 16.11 the perturbed operator L is similar to a self-adjoint operator and its 
disjoint spectral components a/ Cl, i = 0, 1, satisfy the inclusions ( 15.21 ) with ry given by 
(15. 19b - The latter implies that, for i = 0, 1, 



and hence, by ( 15.9 



£/ / d^ 

8i = dist(a,-, a(_,-) > d - r v = - + J — - \\V\ 



\V\\ \\V\\ (\ 21171 

— — < - = tan - arcsin — — - I . (6.6) 



In addition, Theorem 16. II also shows that the spectral subspaces associated with a/ are max- 
imal uniformly definite. Thus all assumptions of Theorem 11.21 (i) (11.21 (ii), respectively) are 
satisfied and the claimed bound (16.31 ) ( (16.51 ). respectively) follows from (11.171 ) ( dl.191 ), respec- 
tively) together with (16.61) . □ 



Remark 6.5. The operators tan©,, i = 0,1, in Theorem 16.41 coincide with the moduli or 
\K*\ of uniformly contractive solutions to the Riccati equations (11.1 lb and (12.4b . respectively 
(see Remark I2T61 and Lemma l2!9l : hence we always have tan©,- < 1. 

If L = A + V satisfies condition (i) in Theorem 16. II i.e. ||V|| <dj%, then the bound (16.3b is 
always less than 1 ; if L satisfies condition (ii) in this theorem and hence ||V|| <d/2, then for the 
bound (16.3b to be less than 1 we need to have d and V such that 

„ „ d ( 2\ d An d 

V < - sin 2arctan- = ^ w - x 0.9060367012. 

11 11 2 V 71 J 24 + 7T 2 2 

Remark 6.6. The bound (16.3b is stronger than the previously known bound 

(\ 7C||V||\ (\ rc||V||\ 
tan©,- < tanh - arctanh = tan - arcsin , i = 0, 1, (6.7) 



2 d J \2 d 

from J6l Theorem 5.8 (i)] and extends it to perturbations V that do not satisfy the condition 
||V|| < dj% required therein. The former is a consequence of the trigonometric inequality 
f tan (j arcsin(2f)) < tanh (^arctanh^?)), t £ (0,1/k). Note that, if ||V|| < d/n, then the 
bound (16.7b may be written equivalently as 

TrllVll 

sin20,-<-V i 5 i = 0,l. (6.8) 
d 



For the particular case (16.4b . the bound (16.5b coincides with the previously known bound 
(11.14b from [6, Theorem 5.8 (ii)] since then the corresponding spectral bounds ry and r (defined 
in (15.13b and (15.19b . respectively) coincide (see Remark I53T ). 
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Remark 6.7. By (15.21) and (15.131) . we also have the estimate 



5 = dist(ffo,Oi) > d-2rv = ^d 2 -4\\V\\ 2 . (6.9) 

Combining estimate (11.211 ) from Theorem ll.3l (i) with inequality ( 16.91 ), we arrive at the bound 

% ||V|| 
tan ®i < = , i = 0,l, 

2 x / t /2_ 4 || V ||2' 

which is worse than (16.3I ). but still better than the estimate (16.71 ) from [|6l Theorem 5.8 (i)]. 
Combining the estimate (11.231 ) from Theorem ll.3l (ii) with the inequality (16.91 ) yields the bound 

llvll 

tanO, < — " " i = 0,l 



3 ± J 



^d 2 -3\\V\\ 2 
which is worse than the estimate (11.141) . 

Remark 6.8. If ||V|| < d/2 and the spectral sets Co and 0\ are bounded and subordinated, 
i.e. conv(ob) nconv(ai) = 0, then combining inequality (16.91 ) with the estimate (11.251 ) from 
Theorem 1 1.3l (iii) results exactly in the a priori sharp norm bound (11.151) from Theorem ll.il 

7. Quantum harmonic oscillator under a ^^-symmetric perturbation 

In this section we apply the results of the previous sections to the A^-dimensional isotropic 
quantum harmonic oscillator under a & ^-symmetric perturbation. 

Let^)=L 2 (]R Ar )for some iVGN. Assuming that the units are chosen such that h = m = (O = 1, 
the Hamiltonian of the isotropic quantum harmonic oscillator is given by 

(Af)(x) = -Uf(x) + hx\ 2 f(x), Dom(A) = (/GW 2 2 (R A ') / dx \x\ 4 \f(x)\ 2 < oo\, (7.1) 
2 2 [ Jrn J 

where A is the Laplacian and W 2 (K ) stands for the Sobolev space of L 2 (R N ) -functions that 
have their second partial derivatives in /^(M^). 

It is well-known that the Hamiltonian A is a self-adjoint operator in /^(M^) and its spectrum 
consists of eigenvalues of the form 

X n =n + N/2, 11 = 0,1,2,..., (7.2) 



whose multiplicities \i n are given by the binomial coefficients (see, e.g., EBB and the references 
therein) 

At«=( Ar + "" 1 ), n = 0,l,2,.... (7.3) 

For n even, the corresponding eigenfunctions f(x) are symmetric with respect to space reflection 
x i — ► — x (i.e. f(—x) = f{x)). For n odd, the eigenfunctions are anti-symmetric (i.e. f(—x) = 
—f{x)). Hence if we partition the spectrum spec(A) = CoU di with 

(7 = {n + N/2 | n = 0,2,4,...}, Oi = {n + N/2 | n= 1,3,5...}, 

then the subspaces 

of symmetric and anti-symmetric functions are the complementary spectral subspaces of A cor- 
responding to the spectral components Oq and G\ , respectively. Obviously, 

d = dist(cTo,CTi) = 1. 
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Let 3* be the parity operator onZ^R^), (^/)(— x) = f(—x), and 2f the (antilinear) operator 
of complex conjugation, (J?f)(x) = f(x), f 6 /^(R^). An operator V on L2(R W ) is called 
^ ^-symmetric if it commutes with the product 2? , i.e. 

@>3-V = V@>3-. (7.5) 

Clearly, the parity operator ^ is a self-adjoint involution on /^(M^) whose spectral subspaces 

RanE^({+l}) =L 2 , even (R JV ), RanE^({-l}) =L 2 , odd (R Ar ) 

coincide with the respective spectral subspaces (17.41) of the Hamiltonian (I7.ll ). 
From now on, let V be the multiplication operator by a function of the form 

V(x)=ib(x), iGK", 

where b £ Lbo(R^) is real-valued and anti-symmetric, i.e. b(x) £ R and b(—x) = —b{x) for a.e. 
x £ R w . Such an operator V is not only ^ ^-symmetric on /^(R^) (see, e.g., lfT2l Section 3]) 
but also /-self-adjoint with respect to the involution J = Moreover, it is anticommuting with 
& which means that such a V is off-diagonal with respect to the decomposition $j = f)o © fj i . 
By |[T2l Theorem 3.2] the spectrum of the perturbed Hamiltonian L = A + V given by 

(Lf)(x) = ~Af(x) + l -\x\ 2 f{x)+ib(x)f(x), Dom(L) = Dom(A), 

remains real (and discrete) whenever ||V|| = ||fc||«> < 1/2. Furthermore, |[T2l Proposition 3.5] 
implies that for ||V|| < 1/2 the closed 0\\y\\ -neighbourhood of the eigenvalue (17.2b of A 
contains exactly \i n (real) eigenvalues X' nk ,k= l,2,...,/x„, of L, counted with multiplicities, 
where \i n is given by ( 17.31 ). Combining |[T2l Theorem 3.2 and Proposition 3.5] with Theorem |5.4| 
ensures that, in fact, the eigenvalues X' n k satisfy the estimates 

\X^-(n+N/2)\<rv, n = 0,l,2,..., k = l,2,...,/x n , 

where 

ry = 1 1 1 1 oo tan arcsin(2||Z?|| 00 )^ < ||6||oo. 

Further, assume that the stronger inequality ||V|| = \\b\\cc < 1/tt holds. In this case it follows 
from (6) Theorem 5.8 (i)] that L is similar to a self-adjoint operator. At the same time, Theo- 
rem [63] implies the following bound on the variation of the spectral subspaces (17.41) : 

tan@ y < ^tan Qai-csin(2||Z7||oc)^ < 1, 7 = 0,1, (7.6) 

where &j = &(Sjj,Sj'j) denotes the operator angle between the subspace Sjj and the spectral 
subspace S)'j of L associated with the spectral subset aj = spec(L) n O rv ((Jj), j = 0, 1. The 
estimate (17.61 ) improves the corresponding bound for the one-dimensional case obtained in (6j 
Section 6] (cf. Remark IdTBI ). 
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